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ABSTRACT 

The  link  center  of  a  simple  polygon  P  is  the  set  of  points  x  inside  P  at 
which  the  maximal  link-distance  from  x  to  any  other  point  in  P  is  minim- 
ized. Here  the  link  distance  between  two  points  x,y  inside  P  is  defined  to 
be  the  smallest  number  of  straight  edges  in  a  polygonal  path  inside  P  con- 
necting X  to  ^.  We  prove  several  geometric  properties  of  the  link  center  and 
present  an  algorithm  that  calculates  this  set  in  time  0(.n^),  where  n  is  the 
number  of  sides  of  P.  We  also  give  an  0(n  log  n)  algorithm  for  finding  a 
point  X  in  an  approximate  link  center,  namely  the  maximal  link  distance 
from  X  to  any  point  in  P  is  at  most  one  more  than  the  value  attained  from 
the  link  center. 


1.   Introduction 

This  paper  concerns  itself  with  analysis  and  calculation  of  the  link  center  of  a  simple 
planar  polygon  P  having  n  sides.  The  notion  of  a  link  distance  between  two  points  x,y 
inside  P  has  been  recently  introduced  in  [Sul];  it  is  defined  as  the  smallest  number  of 
"links"  (i.e.  straight  segments)  in  a  polygonal  path  connecting  x  and  y  within  P,  and  is  a 
useful  metric  for  path  planning  within  P  when  straight  motion  is  easy  to  accomplish  but 
turns  are  expensive. 
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'"  Suri  [Sul]  presents  a  linear-time  algorithm  for  calculating  the  link  distance  between 
any  two  given  points  within  P,  provided  a  triangulation  of  P  is  given  (such  a  triangulation 
can  be  calculated  in  time  0(n  log  n)  by  the  technique  of  [GJPT],  or  in  almost  linear  time 
by  the  recent  algorithm  of  [TV]).  In  fact,  the  technique  in  [Sul]  provides  a  linear-time 
algorithm  for  decomposing  (a  triangulated)  P  into  k -visibility  regions  for  a  fixed  source 
point  u,  where  the  i-th  visibility  region  consists  of  all  points  within  P  whose  link  distance 
from  V  is  exactly  i.  (In  particular,  the  first  visibility  region  is  just  the  set  of  all  points  in  P 
visible  from  u;  see  [AT2],  [EA],  [CG],  [GHLST]  for  analysis  of  visibility  within  a  simple 
polygon.) 

Suri  [Su2]  has  also  considered  the  problem  of  calculating  the  link  diameter  of  P, 
where  this  quantity  is  defined  as  the  maximal  link  distance  between  any  two  points  in  P. 
He  presents  an  algorithm  for  calculating  this  diameter  in  time  0{n  log  n),  using  an 
interesting  divide-and-conquer  approach  on  the  set  of  vertices  of  P.  In  this  paper  we  con- 
sider the  somewhat  more  difficult  problem  of  calculating  the  link  center  of  P,  defined  to  be 
the  set  of  all  points  x  within  P  whose  m'axitnaf  link  distance  to  any  point  of  P  is  the  smal- 
lest possible.  As  will  be  shown  bclo'v,  the  link  center  is  in  general  not  a  singleton; 
nevertheless  it  is  connected  and  has  certain  convexity-like  properties.  From  the  point  of 
view  of  geometric  location  theory,  any  point  in  the  link  center  can  serve  as  a  location  for  a 
mobile  unit  that  has  to  reach  any  point  within  P  so  that  the  maximal  number  of  turns  it 
needs  to  perform  is  minimized;  alternatively,  such  a  point  can  serve  as  a  location  for  a 
transmitter  that  can  broadcast  to  any  point  within  P  (along  a  path  fully  contained  in  P) 
such  that  the  maximal  number  of  relays  necessary  to  reach  any  point  in  P  is  minimized. 

Section  2  of  this  paper  introduces  the  notion  of  geodesic  and  link  geodesic  paths  along 
with  their  associated  metrics  and  establishes  several  useful  results  about  these  concepts. 
Section  3  contains  an  analysis  of  various  properties  of  the  link  center  and  some  additional 
related  results  needed  for  the  calculation  of  the  link  center.  Section  4  presents  an  O(n^) 
time  algorithm  for  calculating  the  link  center.  Section  5  gives  an  0{n  log  n)  time  algo- 
rithm for  finding  a  point  in  an  "approximate"  link  center,  in  the  sense  that  its  maximal 
link  distance  to  any  point  inside  P  is  at  most  one  larger  than  the  smallest  possible  such 
value,  attained  at  points  in  the  actual  link  center. 


2.   Geodesies  and  Link  Geodesies 

Let  P  be  (the  set  of  points  lyin"  in  the  interior  or  on  the  boundary  of)  a  simple 
polygon  in  the  plane.  A  polygonal  path  between  points  x  and  jc'  in  P  will  be  denoted  by 
listing  the  endpoints  of  the  segments  which  constitute  the  path,  for  example: 
X  =  (i=Xo.  •  •  ■  .**=*')■  If  V  and  w  are  two  vertices  of  P  then  b[v,w]  will  refer  to  the 
polygonal  path  obtained  by  traversing  the  boundary  of  P  clockwise  from  i;  to  u;.  Two 
different  metrics  will  be  used  in  what  follows:  the  Euclidean  metric  and  the  link  metric. 
Given  points  x  and  x'  in  a  polygon  P  and  a  continuous  path  x  from  x  to  x'  (in  fact,  it 


suffices  to  consider  only  polygonal  paths),  the  Euclidean  length  of  the  path  is  defined  in  the 
obvious  way  and  is  denoted  by  Ie{x).  Given  a  polygonal  path  x  =  (x=xo,  ■  .  ■  ,xi,=x')  from 
X  to  x',  the  link  length  of  the  path  is  defined  to  be  the  number  k  of  segments  in  the  path 
and  is  denoted  by  l[^(x).  The  Euclidean  distance  cfgix.x')  between  points  x  and  i'  in  a 
polygon  P  is  defined  by 

dgix^')  =  min  {Igix)} 

over  all  continuous  (polygonal)  paths  x  from  x  to  x'.  The  link  distance  d^ix^')  between  x 
and  x'  is  given  by  .     - 

di{x,x  )  =  mm  {//.(x)} 

over  all  polygonal  paths  x  from  x  to  x'.  (It  is  easily  checked  tljat  d^  is  indeed  a  metric  on 
P.)  We  will  call  a  path  x  connecting  points  x  and  x,'  wjthin  P  g?odesic  if  /^(x)  =  dgix^'). 
Similarly,  we  call  x  link-geodesic  if  //,(x)  =  d£,(x,x').  It  is  easy  to  see  that  there  is  a  unique 
geodesic  path  from  x  to  x'  and  that  this  path  is  polygonal;  this  path  will  be  denoted  by 
g(x,x').  On  the  other  hand,  there  may  be  many  link-geodesic  paths  between  two  points.  A 
path  X  from  x  to  x'  is  called  an  optimal  link-geodesic  path  if  x  is  a  link-geodesic  path  and 
Ie(x)  is  minimum  over  all  link-geodesic  paths  from  x  to  x'.  (It  is  easily  checked  that 
optimal  link-geodesic  paths  are  simple.) 

Concerning  these  concepts,  note  that  the  number  of  links  on  a  geodesic  path  between 
two  points  X,  x'  in  P  can  be  considerably  larger  than  the  link  distance  between  x  and  x' 
(see  Fig.  1).  It  also  turns  that  there  can  be  more  than  one  optimal  link-geodesic  paths 
between  two  points.  Moreover,  two  such  paths  may  fail  to  have  a  deformation  that  takes 
one  of  them  to  the  other  through  a  family  of  polygonal  link-geodesic  paths  (see  Fig.  2). 

In  the  following  we  give  a  few  lemmas  about  basic  properties  of  geodesic  and  link- 
geodesic  paths. 

Lemma  A:  Given  line  segments  xy  and  x'y'  in  P  and  optimal  link-geodesic  paths 
X  =  (x=xo,  .  .  .  ,x4=x')  and  y  =  {y=yo,  ■  ■  ■  ,ym-y'\  any  point  17  on  x  can  see  some  point 
r  on  3^  (i.e.  the  segment  qr  is  contained  in  P). 

Proof:  Let  x,x, +  1  be  a  link  of  x;  if  J  crosses  x,Xn.i  then  every  point  on  x,x, +  i  sees  a  point 
on  y.  Assuming  that  y  does  not  cross  x,x,  +  i,  consider  the  chain  ^o.-'^o.J^i.  •  •  •  <Xk<ym-  Form  a 
simple  polygon  Pi  as  follows:  starting  at  x,,  traverse  the  chain  in  each  direction  until  a 
point  on  y  is  reached,  then  connect  these  two  points  using  the  portion  of  y  between  them. 
The  polygon  P;  contains  x,x, +  ,  and  every  interior  point  of  this  segment  must  see  some  ver- 
tex of  Pi  other  than  X;  and  x, +  1.    By  the  optimality  of  x,  no  interior  point  of  x,Xn.i  can  see 
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another  vertex  of  x,  so  the  vertex  of  Pi  seen  by  an  interior  point  of  XiX,  +  i  must  be  a  point 
on  (although  not  necessarily  a  vertex  of)  y.  As  for  the  points  x,  and  x,  +  i  themselves,  if  for 
example  the  other  neighbor  of  x,  along  the  boundary  of  Pj  is  a  point  of  y  we  are  done;  if 
not,  then  either  x,  is  a  reflex  vertex  of  Pi  or  x,  is  a  convex  vertex  of  Pi  but  then  the  seg- 
ment x,_ix,  +  i  is  not  entirely  contained  in  Pi.  In  either  case  x,  must  see  some  vertex  of  Pi 
other  than  its  two  neighbors  and  by  the  optimality  of  x  this  vertex  must  be  a  point  of  y. 
Since  the  same  is  true  for  x,  +  i,  the  proof  is  complete.  Q.E.D. 

Lemma  B:  The  conclusion  of  Lemma  A  holds  if  x  and  y  are  geodesic  paths  (in  the 
Euclidean  metric)  instead  of  optimal  link-geodesic  paths. 


Proof:  Same  as  above.  Q.E.D. 


w 


Lemma  C:  Given  any  points  x,x' ,yiV  and  a  geodesic  (resp.  an  optimal  link-geodesic)  path 
X  =  (x=xo,  .  .  .  ,x^=x')  in  P  and  given  any  point  q  on  x,  one  has 

di(q,y)^mzx  {di{x,y),di{x' ,y)). 


Proof:  Assume  that  x  =  (x  =  xo X4  =  x')  is  a  geodesic  (resp.  an  optimal  link-geodesic) 

path  from  x  to  x'  and  let 

d  =  max  {d/,(x,^),cit(x',y)}. 

We  proceed  by  induction  on  d.  The  case  cf  =  1  is  simple  to  verify.  Indeed,  in  this  case  both 
segments  x^-,  x'^  are  contained  in  P,  and  x  must  be  outward  convex  relative  to  these  two 
segments,  as  is  easily  verified.  In  particular  each  point  on  x  is  visible  from  y  (if  y  coincides 
with  either  x  or  x'  the  argument  is  even  simpler).  Thus  assume  that  ii>l.  Assume  also 
without  loss  of  generality  that  d^ixj)  =  d  and  that  di(x',y)  =  k  where  l^k^d  (the  case 
A  =  0  is  trivial)  and  let  u  =  {y  =  UQ,  .  .  .  ,Ud  =  x)  and  U  =  {y  =  uo,  .  .  .  ,vi,=x')  be  paths  real- 
izing these  distances.  In  addition,  let  uJ  be  a  geodesic  (resp.  an  optimal  link-geodesic)  path 
from  uj_i  to  U;t-i-    By  induction,  every  point  r  on  u;  satisfies 

di{r,y)^ma.x  {dL{ud-i,y),  di{Uk-i^)}  =  d  -  I. 

By  Lemma  B  (resp.  Lemma  A),  each  point  q  on  x  sees  some  r  on  w  and  so 

dc{q,y)^dL(q,r)  +  d^irj)  <  d. 

Q.E.D. 

Lemma  D:  For  any  points  y,y'(.P,  there  exists  a  vertex  y  of  P  such  that  diiy ,y')^di(y ,v). 

Proof:  Let  k  =  di(y,y')  and  let  y  -  (y=yQ,  .  .  .  ,yi,=y')  be  an  optimal  link-geodesic  path 
from  y  to  y'.  Let  w  be  the  point  obtained  by  intersecting  the  boundary  of  P  with  the  ray 
emanating  from  yk-i  which  passes  through  yi,  =  y',  and  let  v  and  u'  be  the  endpoints  of 


the  boundary  edge  thath  contains  w. 

The  segment  connecting  y^-x  and  w  is  certainly  the  optimal  link-geodesic  path  between 
these  two  points.  Applying  Lemma  C  to  this  path  (with  y'  playing  the  role  of  q)  implies 
that  di{y,w)^di{y,y').  Again  applying  Lemma  C,  but  now  to  the  optimal  link-geodesic 
path  formed  by  the  edge  connecting  v  and  v'  (with  w  playing  the  role  of  q),  then  yields 

max{d[^{y,v),di{y,v')}  >  di{y,w)  >  di^{y,y'), 

as  desired.   Q.E.D. 


3.   The  Link  Center 

In  this  section  we  analyze  several  properties  of  the  link  center  that  will  be  useful  for 
the  calculation  of  this  set,  to  be  presented  in  the  following  section. 

The  k-neighborhood  or  k-disk  about  a  point  xiP  (closely  related  to  the  notion  of  i-th 
visibility  regions  mentioned  in  the  introduction)  is  defined  by 

N^(x)  =  {x'iP  I  dL{x^')^k] 

and  the  covering  radius  c(x)  of  x  is  the  smallest  k  such  that  PCNi,{x).  (Note  that  cix)  =  1 
if  and  only  if  P  is  star-shaped  with  respect  to  x.)  The  link  diameter  of  P  is  defined  by 
Dr(P)  =  max  c(x)  and  the  link  radius  is  defined  by  r^(P)  =  min  c(x).  The  link  center  of  P 

X«P  ZiP 

is  defined  by 

Ct(P)  =  {x(P  I  c(x)-=  r^iP)}  . 
See  Figure  3  for  an  illustration  of  these  concepts. 

Using  the  Euclidean  metric  one  can  similarly  define  for  any  non-negative  number  a  the  a- 
neighborhood  about  a  point,  the  covering  radius  of  a  point,  and  the  geodesic  diameter, 
radius  and  center  of  a  polygon.  As  it  turns  out,  these  concepts  for  the  continuous  Euclidean 
metric  are  simpler  to  analyze  and  calculate,  as  demonstrated  e.g.  in  [AT],  [Su2],  [PS]. 

We  will  now  establish  some  properties  of  the  A-disks  Ni,{x)  for  points  xfP  and  of  the 
link  center  of  P. 

Theorem  1:  Let  k^l  and  let  d  be  a  point  of  P. 

The  neighborhood  Ni,(v)  is  a  subpolygon  of  P  all  whose  corners  lie  on  the  boundary  of  P 

and  each  of  whose  edges  is  of  one  of  two  types: 

(a)  a  ^portion  of  a)  side  of  P; 

(b)  a  segment  wz  connecting  a  reflex  vertex  w  of  P  to  a  point  z  on  the  boundary  of  P. 


N2  /         Nj 


-    -- 

laJQD 

■  :■}(.    fc 

-uodis 

■oii-.- 

.    JjClv 

3  Cc  ) 


•••.  \- 


LINK 


V 


/ 


***-.\-i=C  ENTER /y' 


rc^.  ui) 


-6- 


Proof:  (by  induction  on  k) 

For  the  base  case  k  =  l  note  that  Niiv)  is  nothing  but  the  visibility  polygon  of  the  vertex  u, 
which  obviously  has  the  properties  prescribed  in  the  theorem  (see  e.g.  [CG]). 
Now  let  ;fe  >1  and  assume  Ni,-i{u)  has  the  prescribed  properties.  Let  e  be  a  tjrpe  (b)  edge  of 
Nii^iiv).  It  cuts  P  into  two  parts,  one  of  which  does  not  contain  Ni,-i{v).  Call  that  part  P,. 
Let  Vf  be  the  weak  visibility  polygon  of  edge  e  in  polygon  P,,  i.e.  V,  comprises  all  points 
x€P,  for  which  the  straight  segment  xy  is  contained  in  P^  for  some  y^e.  Again,  it  is  well 
known  [CG]  that  V,  is  a  subpolygon  of  P^  that  has  the  structure  prescribed  by  the  theorem 
(with  respect  fo  P,). 

We  claim  that  iV»(u)  =  Ni,.i{v)U[J{V,\e  is  type  (b)  edge  of  iV*_i(i;)}.  This  immediately 
impUes  that  Ni,{u)  has  the  desired  properties,  since  the  boundaries  of  Ni,-i{v)  and  the  visi- 
bility polygons  V,  consist  only  of  type  (a)  and  (b)  edges  and  since  V,  and  Nif-iiu)  intersect 
in  e. 

For  a  proof  of  the  claim  note  that  obviousif"iV*-i(i/}  C  Ni,{v)  and  V,  C  N^(v)  for  each  type 
(b)  edge  of  N't_i(t;).  On  the  other  hand  let  x^Nkiv):  If  diiv^Xk,  then  x^Nk-iiv).  If 
di(.v,x)  =  k,  then  let  x  =  (v=Xo,  .  .  .  ,xi,=x)  be  a  link-geodesic  path  from  i;  to  i.  Clearly, 
Xk-iiNk-iiv)  and  x  not  in  N^t_i(u).  Thus  the  link  x^.ix  has  to  intersect  the  boundary  of 
Nif-iiv),  and,  since  x  lies  entirely  within  P,  this  intersection  has  to  be  at  a  type  (b)  edge  e 
of  Njf^iiv).  Obviously,  xiV^,  which  concludes  the  proof  of  the  claim  and  the  theorem. 
Q.E.D. 

Call  a  subset  SCP  geodesically  convex  if  for  any  points  x,x'iS ,  the  geodesic  path  from 
I  to  x'  is  contained  in  S  (cf.  also  [PS]);  similarly,  call  S  link  convex  if  any  optimal  link- 
geodesic  path  from  x  to  x'  is  contained  in  S.  Note  that  these  notions  of  convexity  are 
preserved  under  intersection.  It  can  be  shown  that  link  convexity  implies  geodesic  convex- 
ity, whereas  the  converse  is  not  true  (see  Fig.  4). 

Lemma  2:  For  any  yiP  and  any  k^O,  Niiiy)  is  both  geodesically  convex  and  link  convex. 

Proof:  This  is  an  immediate  consequence  of  Lemma  C.  Q.E.D. 

Since  the  link  center  of  P  is  just  the  intersection  of  the  sets  N^ix)  over  all  x€P  where 
r  —  ri,(P)  and  since,  by  definition,  the  intersection  of  link  convex  (geodesically  convex)  sets 
is  link  convex  (geodesically  convex),  we  obtain  the  following  result: 

Theorem  3:  lue  link  center  of  a  polygon  P  is  link  convex,  geodesically  convex  and  thus 
connected. 

As  was  mentioned  above,  the  link  center  of  a  polygon  P  consists  of  the  intersection  of 
the  sets  iVr(x)  over  all  xeP.   In  fact,  a  slightly  stronger  result  can  be  proven: 


fu'j     dc  It  r 
r  h '       ■    ■ 
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Theorem  4:  The  link  center  of  a  polygon  P  is  the  intersection  of  the  sets  N^iv)  over  all  con- 
vex vertices  v  of  P,  where  r  =  r£,(P). 

Proof:  The  theorem  clearly  follows  from  Lemma  D  if  we  consider  the  intersection  of  Nr{v) 
for  all  vertices  v  of  P.  However,  Lemma  C  is  easily  seen  to  imply  that  for  a  reflex  vertex  v 
the  neighborhood  N^iv)  contains  the  intersection  of  N,.{vi)  and  Nr{v2),  where  Ui  and  U2  are 
the  two  convex  vertices  at  the  ends  of  the  chain  of  reflex  vertices  of  P  that  contains  v. 
Q.E.D. 

In  the  remainder  of  this  section  we  establish  a  relationship  between  the  link  radius 
and  the  link  diameter  of  a  polygon.   To  accompUsh  this  we  use  the  following  lemmas: 

Lemma  5:  Let  T  be  a  finite  free  tree  and  let  T  be  a  family  of  subtrees  of  T. 
If  rnT'*0  for  every  T,T(:T,th€n  f]{TiT}=fiZ. 

Proof:  See  [G,  Prop.  4.7,  p.92].   Q.E.D.  ^.a  j    .  ..!",r 

Lemma  6:  Let  ^  ^1  and  let  S  be  a  finite  set  of  points  of  polygon  P. 
UNi,(x)nN^{y)*0  for  every  x,:y6S,  then  ni-^*  +  i(^)l*^S}  =6  0. 

Proof:  Let  A  be  an  arbitrary  but  fixed  triangulation  of  P.  The  triangles  of  A  with  the  rela- 
tion "sharing  an  edge"  form  a  graph  that  is  a  tree  T.  For  every  i€S  let  T^  be  the  triangles 
of  A  that  intersect  Ni^iz).  T^  viewed  as  a  graph  must  be  a  subtree  of  T.  The  assumptions 
of  the  lemma  imply  that  T^  and  Ty  have  a  non-empty  intersection  for  every  x,y^S.  Lemma 
5  now  implies  that  P)(T'ilx€S}  *  0,  i.e.  there  must  be  a  triangle  in  A  that  intersects 
Ni,(x)  for  all  x€S.  Clearly  this  triangle  is  contained  in  Ni,+i{x}  for  each  x€S,  which  proves 
our  lemma.  Q.E.D. 

We  are  now  ready  to  prove  that  the  link  radius  of  a  simple  polygon  is  approximately 
half  of  the  link  diameter. 


Theorem  7:  For  any  simple  polygon  P 


Dl<Pi 


2 


^  riCP)  £ 


Dl(P) 


+  1. 


Proof:  Let  ^  =  pi,(P)/2  .    By  the  triangle  inequality  Di{P)<2ri^{P),  i.e.    k^r^^iP).    We 

claim  that  Nii{x}nNi,(y)^0  for  any  pair  of  points  i,^€P,  which  by  virtue  of  Lemma  6  and 
Theorem  4  implies  the  desired  r[^{P)<k  +1.  For  a  proof  of  the  claim  let  r  and  y  be  two 
points  in  P.    Let  w  =  {x  —  wq,  .  .  .  ,Ws=y)  he  a  link-geodesic  path  connecting  x  and  y  and  let 
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i=\s/2].  Since  s  cannot  exceed  the  diameter  D/,(P)  we  have  i^k.  Consider  u;,:  Obviously 
(fi,(i,u;,)<i  and  di{w„y)^s -i^i^k.  Thus  the  neighborhoods  N,{x)  and  N,{y}  intersect 
and  so  do  N^lx)  and  Ni,(y},  which  proves  the  claim  and  the  theorem.  Q.E.D. 

It  is  possible  to  strengthen  the  result  of  Theorem  7  somewhat.  For  polygons  with  odd 
link-diameter  DiiP)  we  can  show  that  the  link  radius  r/,{P)  is  always  equal  to  pi,(P)/2  . 
The  proof  of  this  fact  is  very  similar  to  the  proof  given  above  but  it  involves  more  geometry. 

At  first  we  give  an  analogue  of  Lemma  6.  It  is  a  topological  version  of  Helly's 
theorem  for  convex  sets,  first  shown  by  Helly  himself  [H]  and  then  in  a  slightly  stronger 
form  by  Molnar  [M].  For  our  purposes  a  cell  in  the  plane  is  a  simply  connected  compact 
subset  of  the  plane. 

Lemma  8:  Let  C  be  a  set  of  cells  in  the  plane.  ■  *  ^^'^  '----'^' 

If  CnC  is  a  cell  for  every  C,C'€C  and  CnC'nC'^^  0  for  G,CsCf 'eC,  theV  C\{C^C}^0. 

With  this  lemma  we  can  prove  a  slightly  stronger  version  of  Lemma  6. 

Lemma  9:  Let  k>l  and  let  S  be  a  finite  set  of  points  of  polygon  P. 

U Nk(x}nNk  +  i{y)^0  foT  every  x,y(:S ,  then  f]{Ni,  +  i(x)\x^S}  *  0. 

Proof:  First  note  that  clearly  all  neighborhoods  iV^(x)  are  cells.  Any  two  such  neighbor- 
hoods, if  they  have  non-empty  intersection,  intersect  in  a  cell  (because  they  are  link  con- 
vex). Thus  by  Lemma  8  it  suffices  to  show  that  the  assumptions  of  this  lemma  imply  that 
iVi  +  i(x)n^*  +  iO')n^i  +  ,(2)*0forx,y,2€S. 

Let  x,y^  be  three  points  in  S  and  consider  the  neighborhoods  ^^  +  i(x),  iV^  +  i(y),  and 
N^iz).  U  Ni,^i(x)r\Ni,^i{y)nN^{z)^0,  then  obviously  also  Ni,^iix)nN^  +  i(y)nN,,^i(z)*0, 
as  desired. 

So  assume  that  Ni,^.i(.x)r\N)i  +  i{y)nNi,{z)  =  0.  The  assumptions  of  this  lemma  imply  that 
^*  +  i(^)'^-^*  +  iCy)*25.  Because  of  link  convexity  (Lemma  2)  this  intersetion  must  be  sim- 
ply connected.  Thus  Ni,  +  i{x}UNii  +  i(y)  forms  a  subpolygon  of  P  that  is  simple,  i.e.  it  has  no 
holes.  By  the  assumptions  of  this  lemma  Ni,{z)  intersects  both  7V^  +  i(x)  and  Ni,  +  i{y).  Since 
we  assume  that  Ni,  +  i{x)nNi,  +  i(y)nNi,{z)  =  0,  it  follows  that  iV^(z)  intersects 
Nit  +  i(x)UNi,  +  i(y)  in  at  least  two  different  connected  components.  Thus  the  complement  of 
N|,  +  l{x)UN|,^.l(y)UN|,{z)  has  a  bounded  connected  component  (i.e.  the  polygon  formed  by 
this  union  has  a  hole).    Call  it  X. 

Consider  the  boundary  of  X.  Since  X  lies  completely  inside  the  polygon  P,  the  struc- 
ture Theorem  1  implies  that  each  of  Ni,  +  i(x),  Ni,  +  i{y),  and  Ni,(z)  can  contribute  at  most  one 
edge  to  the  boundary  of  X.  It  follows  that  X  must  be  a  triangle  and  the  edge  contributed 
by  Ni,(z)  can  "see"  all  of  X,  in  particular  also  its  opposite  corner  (which  is  in 
Ni,Jri(x)C\Ni,^.i{y)).       Thus      this      corner      is      in      Ni,^.i{z)      and      it      follows      that 


iV4  +  i(i)niy^4  +  i(^)niVi  +  i(z)*0,  as  desired.  Q.E.D. 

With  this  lemma  the  following  theorem  can  be  proven  in  the  same  manner  as  Theorem  7. 

Theorem  10:  For  a  simple  polygon  P  with  odd  link  diameter  Di(P)  the  link  radius  is  given 
by 


riCP)  = 


DciP) 


4.  Calculating  the  Link  Center 

The  results  of  the  preceding  section  lead  to  a  rather  simple  approach  to  the  calcula- 

T;'*    f    SCI'.  -. 

tion  of  the  link  center  of  a  given  simple  polygon  P  with  n  sides,  which  can  be  implemented 
to  run  in  0(nh  time.  Specifically,  we  first  calculate  the  link  diameter  of  P  in  time 
0{n  log  n)  by  the  technique  in  [Sul].  Theorems  7  and  10  then  give  at  most  two  possible 
values  for  the  link-radius  r  -  ri^(P):  It  can  be  I)i,(P)/2  or,  in  the  case  of  even  link  diame- 
ter, it  can  also  be  fD/.(P)/2l  +  1. 

By  Theorem  4  we  need  to  calculate  the  intersection  (~]  Nriv)  over  all  convex  vertices 
y  of  P.  We  shall  try  this  for  the  smaller  of  the  two  candidate  values  for  r  first.  If  the 
intersection  is  non-empty,  then  it  is  the  link-center  of  P  If  it  is  found  to  be  empty,  then 
Theorem  7  guarantees  that  the  intersection  pliV^iu)  will  be  non-empty  for  the  larger  can- 
didate value  of  r  and  this  intersection  will  constitute  the  link-center.  Thus  it  suffices  to 
present  an  O(n^)  algorithm  for  the  construction  of  P)iV/,([;)  over  all  convex  vertices  v  of  P, 
where  k  is  any  fixed  integer  between  1  and  n. 

Theorem  1  implies  that  each  Ni,{v)  is  a  subpolygon  of  P  each  of  whose  sides  is  either 
(a  portion  of)  an  edge  of  P  or  a  segment  e  connecting  a  reflex  vertex  u;  of  P  to  the  boundary 
of  P.  Denote  such  a  segment  by  e^iw).  Clearly,  each  reflex  vertex  u;  of  P  can  have  at  most 
one  segment  e^iw)  emanating  from  it  for  each  (convex)  vertex  u  of  P. 

We  therefore  proceed  as  follows.  We  first  triangulate  P  in  time  Oin  log  n)  as  in 
[GJPT].  For  each  convex  vertex  u  of  P  we  calculate  Ni,(v)  in  0(n)  time,  using  Suri's  algo- 
rithm [Sul],  and  for  each  reflex  vertex  w  we  collect  the  (at  most  n)  segments  e^(w).  We 
can  then  calculate  in  Oin)  time  a  subpolygon  Piw)  of  P  whose  boundary  consists  of  two 
segments  ej,  e2  emanating  from  w,  each  of  which  is  either  one  of  the  segments  e^Xw)  or  a 
side  of  P  incident  to  w,  and  of  a  portion  b[qi,q2]  of  the  boundary  of  P  between  the  two 
other  endpoints  qi,  qi  of  ej,  e2  respectively.  Note  that  the  preceding  discussion  implies 
that  the  link  center  is  the  intersection  (^  P{w)  over  all  reflex  vertices  w  of  P.    (Note  also 

w 

that  if  P  has  no  reflex  vertices  then  the  link  center  of  P  is  clearly  P  itself.) 
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We  are   thus   left   with   the   problem   of  calculating    f)  P{w).     Note   that  the   total 

w 

number  of  edges  in  the  boundaries  of  the  subpolygons  P(m;)  is  only  0(n),  so  that  the  inter- 
section of  these  polygons  can  be  calculated  trivially  in  time  O(ra^)  as  follows.  We  calculate 
the  Oin'^)  intersection  points  between  edges  of  (distinct)  P(w)'s.   For  each  edge  e  of  such  a 


polygon  we  have  Oin)  intersection  points  of  this  kind  along  e,  so  that  e  fl 


n  P('^) 


is 


equal  to  the  intersection  of  half-segments  of  e,  each  delimited  at  one  of  these  intersections. 
The  intersection  of  these  half-segments  is  trivial  to  calculate  in  0(n)  time.  Thus  all  the 
edges  along  the  boundary  of  the  link  center  can  be  calculated  in  overall  O(n^)  time.  Since 
there  are  only  0(n)  such  edges,  the  final  step  of  finding  the  circular  order  of  these  edges 
along  the  boundary  of  the  link  center  is  now  easy  to  accomplish  in  subquadratic  time.  We 
have  thus  shown 

Theorem  11:  The  link  center  of  a  simple  polygon  P  having  n  edges  can  be  calculated  in 
Oin^)  time. 


5.   Finding  an  approximate  link  center 

The  result  of  the  preceding  section  immediately  raises  the  question  whether  the  link 
center  of  a  simple  polygon  P  with  n  sides  can  be  calculated  in  subquadratic  time.  This 
problem  is  also  motivated  by  the  fact  (used  above)  that  the  link  diameter  (and  an  approxi- 
mation of  the  link  radius)  can  be  calculated  in  0{n  log  n)  [Su2],  and  that  the  (unique)  geo- 
desic center  of  P  can  be  calculated  in  time  0{n  log^n)  [PS].  Even  if  all  we  seek  is  just  one 
point  in  the  link  center,  we  still  do  not  know  how  to  find  it  in  subquadratic  time.  However, 
as  we  show  in  this  section,  if  we  are  willing  to  compromise,  and  are  satisfied  with  finding 
an  approximate  link  center,  by  which  we  mean  a  point  whose  covering  radius  exceeds  the 
link  radius  at  most  by  one,  then  we  can  achieve  this  goal  in  time  Oin  log  n). 

Let  A  be  an  arbitrary  triangulation  of  P  and  let  r  =  r/,(P)  be  the  link  radius  of  P.  The 
following  observation  is  almost  too  trivial  to  state  (see  also  the  proof  of  Lemma  6): 

Lemma  12:  There  is  a  triangle  t  in  the  triangulation  A  such  that  every  point  in  t  is  an 
approximate  link  center  of  P. 

In  other  words  every  xdt  has  covering  radius  C£,(x)  <  r  +  l.  Let  us  call  such  a  triangle 
a  center  triangle.  That  a  center  triangle  must  exist  is  clear:  Take  any  triangle  that  con- 
tains a  point  of  the  true  link  center  of  P. 

How  can  we  find  a  center  triangle  quickly?  Some  sort  of  "binary  search"  on  the  tree 
formed  by  the  triangles  of  A  seems  to  be  called  for.    Recall  that,  as  pointed  out  already  in 


11  - 


the  proof  of  Lemma  6,  the  triangles  of  A  form  a  graph  T  that  is  a  tree.  In  this  graph  the 
n-2  nodes  are  the  triangles  of  A  and  two  nodes  are  adjacent  iff  the  corresponding  triangles 
share  an  edge.  Obviously  no  node  in  T  has  degree  exceeding  three.  It  is  well  known  that, 
using  only  linear  time,  one  can  always  find  in  such  a  tree  an  edge  whose  removal  would 
leave  two  subtrees,  each  rintaining  no  more  than  2/3  of  the  nodes  of  the  original  tree.  If 
we  could  also  determine  in  0(n)  time  which  of  the  two  subtrees  contains  a  node  correspond- 
ing to  a  center  triangle  in  A,  then  recursive  application  of  this  tree  splitting  idea  would 
yield  an  0{n  log  n)  algorithm  for  finding  such  a  center  triangle.  (See  [Ch]  and  [GHLST] 
for  other  algorithms  based  on  this  idea.) 

Thus  we  want  to  solve  the  following  problem  in  linear  time:  Given  an  edge  e  of  the 
triangulation  A  that  cuts  the  polygon  P  into  two  polygons  Pi  and  P2,  which  of  Pi  or  P2 
contains  a  triangle  of  A  that  is  a  center  triangle  of  P? 

At  first  a  few  more  definitions:  In  the  usual  way  we  can  generalize  the  distance  func- 
tion di  so  that  one  of  the  arguments  is  a  subset  of  P: 

di(A,x)  =  min{di(y,x)\yiA} 

The  k-neighborhood  of  a  set  ACP  in  P  is  then  defined  as 

iV*(A,P)  =  {x6P|dt(A,x)<fe}, 
and  the  covering  radius  of  A  is  defined  as 

c{A  ,P)  =  mini/:  |  PCN*(A  ,P)}. 
With  these  definitions  we  can  state  the  following  characterization: 

Lemma  13:  Let  e  be  a  diagonal  of  a  triangulation  A  of  a  simple  polygon  P  that  cuts  P  into 
two  simple  polygons  Pi  and  P2.   Let  Ci  =  c(e,Pi)  and  C2  =  c(e,P2). 

(i)       If  Ci=C2,  then  both  triangles  of  A  that  contain  e  are  center  triangles  for  P. 

(ii)     If  Ci<C2,  then  P2  contains  a  triangle  of  A  that  is  a  center  triangle  for  P. 

Proof:  Let  r  =  r^(P)  be  the  link  radius  of  the  polygon  P. 

(i)  Let  2  be  a  point  of  the  link  center  of  P.  Without  loss  of  generality  assume  z  lies  in  P^ 
Let  p  be  a  point  in  P2  for  which  di(p,e)  =  C2.  Since  a  link  geodesic  path  connecting  p  and  z 
must  cross  the  segment  e,  we  have  Ci=C2^d[^(p,z)^r.  From  this  it  follows  immediately 
that  both  triangles  of  A  that  contain  e  must  be  center  triangles  of  P  since  any  point  q€P 
can  be  connected  to  any  point  x  in  these  two  triangles  by  a  link  path  with  at  most  r  +  l  seg- 
ments (at  most  r  segme  Us  from  q  to  some  point  on  e  and  one  more  segment  to  x). 

(ii)  If  P2  contains  a  point  of  the  true  link  center  of  P,  then  there  is  nothing  to  prove  since 
the  triangle  of  A  that  contains  that  point  will  be  a  center  triangle  of  P  contained  in  P2.  So 
assume  Pi  contains  a  point  2  of  the  link  center  of  P.  Using  exactly  the  same  argument  as 
in  (i)  we  get  the  inequalities  Ci<C2^'",  and  as  in  (i)  we  can  conclude  that  the  triangle  of  A 
that  contains  the  edge  e  and  is  contained  in  P2  is  a  center  triangle.  Q.E.D. 
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The  only  question  that  remains  is  whether  the  covering  radii  c(e,Pi)  and  c(e,P2)  can 
be  computed  in  linear  time.  Using  the  techniques  presented  in  [Sul]  this  is  indeed  possible 
if  a  triangulation  A  of  P  is  available.   Thus  we  summarize  the  result  of  this  section: 

Theorem  14:  For  a  simple  polygon  P  with  n  vertices  it  is  possible  to  compute  in  Oin  log  n) 
time  a  triangle  contained  in  P  all  of  whose  points  are  approximate  link  centers  of  P,  i.e. 
they  have  a  covering  radius  that  exceeds  the  link  radius  of  P  at  most  by  one. 

(Note:  The  results  of  this  section  provide  an  alternative  method  for  obtaining  an  approxi- 
mation for  the  link  radius  with  error  at  most  one,  instead  of  calculating  the  link  diameter 
and  using  Theorem  7  or  Theorem  10.) 

References  .       _- 

[AT]  Asano,  A.  and  Toussaint,  G.T.  (1985).  Computing  the  geodesic  center  of  a  simple 
polygon,  Tech,  Report  SOCS-85.32,  McGill  University,  see  also  Proc.  US  -  Japan 
Seminar  on  Computer  Science,  Japan,  June  1986. 

[AT2]  Avis,  D.  and  Toussaint,  G.T.  (1981).  An  optimal  algorithm  for  determining  the 
visibility  of  a  polygon  from  an  edge,  IEEE  Trans.  Computers  C-30,  pp.  910-914. 

[Ch]  Chazelle,  B.  (1982).  A  theorem  on  polygon  cutting  with  applications,  Proc.  23rd 
IEEE  Symp.  on  Foundations  of  Computer  Science,  pp.  339-349. 

[CG]  Chazelle,  B.  and  Guibas,  L.  (1985).  Visibility  and  intersection  problems  in  plane 
geometry,  Proc.  ACM  Symposium  on  Computational  Geometry,  pp.  135-146. 

[EA]  El  Gindy,  H.  and  Avis,  D.  (1981).  A  linear  algorithm  for  computing  the  visibility 
polygon  from  a  point,  J.  Algorithms  2,  pp.  186-197. 

[GJPT]  Garey,  M.R.,  Johnson,  D.S.,  Preparata,  F.P.  and  Tarjan,  R.E.  (1978).  Triangulat- 
ing a  simple  polygon,  Inf.  Proc.  Letters  7,  pp.  175-180. 

[G]  Golumbic,  M.C.,  Algorithmic  Graph  Theory  and  Perfect  Graphs.    Academic  Press 

(1980). 

[H]  Helly,  E.  (1930).    Uber  Systeme  abgeschlossener  Mengen  mit  gemeinschaftlichen 

Punkten,  Monatsh.  f.  Math.  u.  Physik  37,  pp.  281-302. 

[GHLSTlGuibas,  L.,  Hershberger,  J.,  Leven,  D.,  Sharir,  M.,  and  Tarjan,  R.E.  (1986).  Linear 
time  algorithms  for  visibility  and  shortest  path  problems  inside  a  triangulated  sim- 
ple polygon,  Tech.  Report  218,  Comp.  Sci.  Dept.,  Courant  Institute,  April  1986;  (to 
appear  in  Algorithmica). 

[LP]  Lee,  D.T.  and  Preparata,  F.P.  (1984).  Euclidean  shortest  paths  in  the  presence  of 
rectilinear  barriers,  Networks  14  (3),  pp.  393-410. 


13 


[M]  Molnar,  J.  (1957).  Uber  den  zweidimensionalen  topologischen  Satz  von  Helly,  Mat. 
Lapok  8,  pp.  108-114.  (Hungarian  with  German  and  Russian  summaries) 

[PS]  Pollack,  R.  and  Sharir,  M.  (1986).  Computing  the  geodesic  center  of  a  simple 
polygon,  Tech.  Report  231,  Comp.  Sci.  Dept.,  Courant  Institute,  July  1986;  (to 
appear  in  Discrete  and  Computational  Geometry). 

[Sul]  Suri,  S.  (1986).  A  polygon  partitioning  technique  for  link  distance  problems. 
Manuscript,  Dept.  of  Comp.  Science,  Johns  Hopkins  University,  November  1986. 

[Su2]  Suri,  S.  (1986).  Computing  all  geodesic  furthest  neighbors  of  a  simple  polygon. 
Manuscript,  Dept.  of  Comp.  Science,  Johns  Hopkins  University,  November  1986; 
(to  appear  in  Proc.  of  the  3rd  ACM  Conference  on  Computational  Geometry). 

[TV]  Tarjan,  R.E.  and  Van  Wyk,  C.J.  (1986).  An  0{n  log  log  n)  time  algorithm  for  tri- 
angulating simple  polygons,  preprint  (submitted  to  SIAM  J.  Computing). 


'^U  COMPscr    n- 


